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PREFACE 


This document contains the new Unified Conic formulation of R. H. Battin, 
and is intended to replace the Space Shuttle GN&C Equation Document No. 3 by 
W. M. Robertson. 
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FOREWORD 


This document is one of a series of candidates for inclusion in a future 
revision of JSC -04217, "Space Shuttle Guidance, Navigation and Control Design 
Equations." The enclosed has been prepared under NAS9-10268, Task No. 15-A, 
"GN&C Flight Equation Specification Support", and applies to function 1 of the 
Orbital Coast Navigation Module (ON2) as defined in MSC-03690, Rev. D, Space 
Shuttle Orbiter Guidance, Navigation and Control Software Functional Require- 
ments", dated January 1973, 



Division Leader, Guidance Analysis 
NASA Programs Department 
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NOMENCLATURE 


a 

A 

^4 

B 

D 


Semi- major axis 

Value of ^ U2(x; or ) function 

Intermediate variables 

Reciprocal of normalized final position magnitude (b^ = ^0^^^ 
Value of 2 (Ug {•^; a; ) + Uq ® ^ ^ function 

*Value of Uj (^ ; or ) function 


®flag Error flag 


E * Value of Uq J o' ) function 

f True anomaly 

fg Coefficient indicating sign of transfer 

f Coefficient to convert to normalized velocity ( v = f y ) 

fg Inverse of f " ^2 — ^ 

fg Coefficient to convert to normalized transfer time interval 

(At = fg At) 

Inverse of f^ (At = f ^ A? ) 

fg Coefficient to convert to normalized eccentric anomaly 

(x = fg X ) 


These capital letters do not indicate matrices. 
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Inverse of fg (x = fg x ) 


^7 

F 


F,,F, 


G,G 

G* 


Intermediate variable 


State extrapolation coefficient 


State extrapolation coefficients (F^ 


"state extrapolation coefficients 


State extrapolation coefficient 


■M 


Ff) 



(G=/^ G ) 


Iteration counter 


max 




3 

k 


Maximum permissible number of iterations 


Unit vector in r ^ direction 


Counter for calculation of the continued fraction Q 


Counter for the summation of the series 


max 

P 

P 

Q 

r, r 
r(t) 


Maximum value for k 


Normalized semi-latus rectum (p = Tq p ) 

C 

Normalized orbital period ( P = Pgi^cular ^ ^ 


Continued fraction function 


Final position magnitude and its normalized analog (r = r^ r ) 


Inertial position vector corresponding to time t 


Magnitude of Tq 


These capital letters do not indicate matrices. 
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Inertial position vector corresponding to time t^ 


max 


Maximum scalar representable on the computer 


®min 


Minimum positive scalar representable on the computer 


u. 


Value of ) function 


u 


lold 


Previous value of u- 


U^(x;a) 


Universal function of Battin 


v(t) 


Inertial velocity vector corresponding to time t 


Zo' “0 Inertial velocity vector, corresponding to initial time t^, 

and its normalized analog (vq = Vp ) 


w 


Value of Uq ; a ) function 


x,x 


Eccentric anomaly difference and its normalized analog 
(x = X , independent variable used in the Kepler 

iteration scheme) 


arg 


Independent variable used as the input argument to the 
Transfer Time Interval Routine 


Value of X to which the Kepler iteration scheme converged 


Previous value of 


‘guess 


Guess of X 


‘last 


Previous value of x 


This capital letter does not indicate a matrix. 
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*max’ ^max 


^min* ^min 


‘old 




z, z 
a , a 

At. aT 
^^arg 

At' 


At 


error 


At 


guess 


At 


last 


AV 

max 


At’ , 
mm 


Independent variable upper bound, and its normalized analog 

Independent variable lower bound, and its normalized analog 

Value of X from the previous iteration 

guess 

Value of X corresponding to one period of the orbit 
Intermediate variables 

Values of Uj , Ui(f- ; a ) functions (z = z ) 

Reciprocal of the semi‘-major axis and its normalized analog 

(a = Tq O' ) 

V 

Transfer time interval and its normalized analog (At “ ^circular 

Dependent variable used as the output argument from the 
Transfer Time Interval Routine 

Value of the transfer time interval calculated in the Kepler Routine 
Previous value of At^ 

Error between ^ and At during any particular iteration 

guess ^ 

Transfer time interval corresponding to the most recent independent 
variable during the Kepler iteration 

Previous value of ^t 

Dependent variable boundary distance on the maximum side used 
in the Secant Iterator 

Dependent variable boundary distance on the minimum side used 
in the Secant Iterator 
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A't, 


max' 

'max 


nun 


aX 


nun 

^^old 


Au 


Ax 


Dependent variable upper bound and its normalized analog 

Dependent variable lower bound and its normalized analog 

Value of ^ from the previous iteration 

guess 

Increment in 
Increment in 

guess 

Tolerance defining near parabolic orbits used in boundary selection 

Primary convergence criterion: error in normalized transfer 
time interval 


Secondary convergence criterion: minimum permissible dif- 
ference between two transfer time intervals used in the Secant 
Iterator 

Tolerance defining area where the computer cannot handle the 
continued fraction Q calculation 

0 Transfer angle (true anomaly increment) 

Gravitational parameter of the earth 

aflj'o'o Conic parameter and its normalized analog ( cFq = ) 


X 
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INTRODUCTION 


The Conic State Extrapolation Routine provides the capabil- 
ity to conically extrapolate any spacecraft inertial state vector either 
backwards or forwards as a function of time or as a function of 
transfer angle. It is merely the coded form of two versions of the 
solution of the two-body differential equations of motion of 
the spacecraft center of mass. Because of its relatively fast compu- 
tation speed and moderate accuracy, it serves as a preliminary 
navigation tool and as a method of obtaining quick solutions for tar- 
geting and guidance functions. More accurate (but slower) results 
are provided by the Precision State Extrapolation Routine. 
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2. FUNCTIONAL FLOW DIAGRAMS 

The Conic State Extrapolation Routine basically consists of two parts - one 
for extrapolating in transfer time interval, and one for extrapolating in transfer 
angle. Some portions of the formulation are common to the two parts and may be 
arranged as subroutines on a computer. Similarly, there are several areas of 
overlap with the Conic Required Velocity Routine (Ref. 5), and these may be treated 
likewise. 

2. 1 Conic State Extrapolation as a Function of 
Transfer Time Interval - Kepler Routine 

This routine involves a single loop iterative procedure, and hence is orga- 
nized into three sections - initialization, iteration, and final computations. This 

scheme is shown in Figure 1, where the variable independent 

guess 

variable used in the iteration procedure. For a given initial state, the variable 

"x " measures the amount of transfer along the extrapolated trajectory, 
guess 

Hence, the iteration involves adjusting " until the transfer time interval 

guess 

calculated from it agrees, to within some tolerance, with the specified transfer 
time interval. Then the extrapolated state vector is calculated from several in- 
termediate variables already calculated in the transfer time interval computation. 

2. 2 Conic State Extrapolation as a Function of 
Transfer Angle - Theta Routine 

This routine involves a direct calculation (i. e. , does not have an iteration 
scheme), as shown in Figure 2. The extrapolated state vector is again calculated 
from the same intermediate variables used in the Kepler Routine. However, this 
time these variables are computed directly in terms of the initial state and the 
transfer angle. 
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ENTER 



EXIT 


Figure 1. Kepler Routine, Functional Flow Diagram 
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ENTER 



EXIT 


Figure 2. Theta Routine, Functional Flow Diagram 
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INPUT AND OUTPUT VARIABLES 


The Conic State Extrapolation Routine has only one universal constant: the 
gravitational parameter of the earth. Its principal input variables are the inertial 
state vector which is to be extrapolated and the transfer time interval or transfer 
angle through which the extrapolation is to be made. Several optional input vari- 
ables may be supplied in the transfer time case in order to speed the computations. 
The principal output variable of both cases is the extrapolated inertial state vector. 

3# 1 Conic State Extrapolation as a Function of 
Transfer Time Interval - Kepler Routine 

Input Variables 

(Fo^Vq) Inertial state vector which is to be extrapolated 
(corresponds to time t^ ). 

At Transfer time interval through which the extra- 

polation is to be made . 

X Starting guess of the independent variable (used to 

speed convergence). If no guess is available, 
set X = 0, and the routine will generate its own 
guess. 

At’ Value of dependent variable, the transfer time inter- 

val, to which the previous call to Kepler had converged. 

Value of the independent variable, universal eccentric 
anomaly difference, to which the previous call to Kepler 
had converged. 


Output Variables 

(r, v) Extrapolated inertial state vector (corresponds to time t). 

At Converged value of the dependent variable, the transfer 

time interval, from the Kepler iteration scheme (should 
agree closely with At), 
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X 


c 


®flag 


Converged value of the independent variable from the 
Kepler iteration scheme* 


Error flag 


! 0 no error 

1 solution failed to converge 
2 series failed to converge 
3 both 1 and 2 occured 


3* 2 Conic State Extrapolation as a Function of 
Transfer Angle - Theta Routine 


(r q,Vq) Inertial state vector which is to be extrapolated* 

0 Transfer angle through which the extrapolation is 

to be made* 


(r , v) Extrapolated inertial state vector • 


Transfer Time Interval corresponding to the conic 
c 

extrapolation through the transfer angle 6. 


®flag 


Error flag = 


0 

1 

2 

3 

4 


no error 

orbit too nearly rectilinear 

multi -revolution input for 
hyperbolic trajectory 

transfer past asymptote on 
a hyperbolic trajectory 

transfer required closure 
through infinity 
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4. 


DESCRIPTION OF EQUATIONS 


The universal formulation of Battin in terms of the universal eccentric anomaly 
difference is used. It should be noted that the formulation given in the detailed flow 
diagrams of Section 5 is a normalized version of the formulation presented here. 

It is felt that the equation descriptions of this section are more clear when ex- 
pressed in standard variables rather than in normalized variables. 

The universal eccentric anomaly, usually denoted by x, is defined by the 
relations 

I ^/T (E - Eq) ellipse 

^/p“(tanf/2 - tanfQ/2 ) parabola 
(H - Hq ) hyperbola 


Here, a is the semi-major axis, E and H are the eccentric anomaly and its 
hyperbolic analog, p is the semi-latus rectum and f the true anomaly. Using this 
definition, the generalized form of Kepler^s equation expresses the transfer time 
interval as 


where 

and 

and 


UflCx ; a) = 


Q Uj(x; a) + a 

1 2 

/ or X 

' n! " (n+2)! 


0 


U2 (x; a) 


2 ^ 
(n+4) ! 



+ Ug(x5 a)] 

- ... ) 


**0' '^0 




Furthermore, the expressions for the extrapolated position and velocity vectors 
(r, V ) in terms of the initial position and velocity vectors (r q , v^ ) are also 
given in terms of x by the following standard formulae 


r(t) 

v(t) 


r = Tq Uq(x; O') + O’ q Uj^(x; a) + U2(x; a) 


= [1 -7— Ugix; ff)] r„ + [ 

0 


-0 


r^j Uj^(x; O') + Oq UjCx; a) 


= [ Uj(x; O') ] Tq + ( 1 Ugix; Of) ] Vj 


^0 
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However, in this document a slightly modified, but equivalent, formulation 
is used, namely 

U (— • a) 

, . . j2(rQ Uq(j ; a) + CTq U^(y; »))+ y U2(x; a)Q] 

T - Tq + (l-or Tq) UgCx; a-) + 2 UQ<'y; a) Uj (y ; o) 

£(t)= FTq+Gvq 

v(t)= F^To+G^Vq 

Here, the function Q is defined in Section 4. 1, but the other variables are as 
follows: 

F = 1 U«(x; a ) 


G = 


2 Tq Uq (y; a) (y; a) + (x; a) 

“ If 


f't ■ -Irf 


= 1 UgCx; a) 


The family of functions or) used in these equations can be related to 

some elementary functions by noting that 


Uj(x; a) = 


\ cos X 

a>0 

|cosh X 

a <0 

/ sin ^ 

a >0 

1 


fsinh J-a x 

a <0 

1 



Then, from these the higher order functions follow directly from the useful identity 


Uj^(x; a) + a " ITT 
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4. 1 The Q Function 


A very useful function employed in the evaluation of the universal formulation 
of Battin is the Q function which may be written 


<x ; or) 
— 

<T ‘ 


It can be shown that Q is a hypergeometric function of the variable 

q = O' a ) 

From this, Battin has developed a continued fraction formulation for evaluating Q 
as a function of the variable 

w = 


to improve the convergence* The formulation is conveniently expressed by 


O = ~ — f 1 +— L 

^ ^2 f 3 w(w+l) 


(1 -T^) ] 


where 


B 


w- 1 
w+1 


n-1 




Evaluation of this expression is started by setting = 0 for a sufficiently large ‘ 
value of M, and calculating lower *s recursively until is reached* This 
means that the number of levels M of the continued fraction necessary to achieve 
some specified accuracy must be precomputed as a function of w. 

4* 2 Conic State Extrapolation as a Function of Transfer 

Time Interval - Kepler Routine 

Since the transfer time interval At is given, it is desired to find the x 
corresponding to it in the generalized Kepler equation* Then the extrapolated state 
vector {r, v) expression can be evaluated based on that value of x . Unfortunately, 
the generalized Kepler equation expresses At as a transcendental function of x , 
and no power series inversion of the equation is known which has good convergence 
properties foi* all orbits* Therefore, it is necessary to solve the equation itera- 
tively for X* 
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The actual evaluation of the Kepler equation is conveniently expressed in terms 
of the variable 


z = 2 Uj if;a) 


From this it can be shown that 

Oo<7‘“> 

Dj<f :. ) 

! » > 

Ug (x; a) 



w z 




These quantities are needed to evaluate the transfer time interval corresponding to 
X during each iteration. Then, once convergence is achieved the latest values of 
these quantities are used to calcxilate the extrapolated state vectors. 


The iteration scheme used to solve for x is a modified secant method (linear 
inverse interpolation/extrapolation). It merely finds a new guess for the independent 
variable which attempts to adjust the dependent variable At to the desired 

value. This adjustment is based on a linear interpolation/ extrapolation of the last 
guess of the independent variable . The scheme requires a minimum point 

and a maximum point on the At vs. x curve between which the solution is known 
to be. These two points define two possible guesses for the new x, namely 


^new 

X 

new 


^old 

^old 


\**old - 
\"old - -i* 


min^ 


max 


<^old ■ ^min> 
^^old ^ ^max^ 


one of which is an interpolation and the other an extrapolation. The modified scheme 
preferentially selects the extrapolated guess from the two possibilities* However, 
if this happens to extrapolate outside the upper or lower bounds, then the interpolated 
guess is used to ensure a solution. It should be noted that these bounds are contin- 
ually reset during the iteration as more and more values of x are determined 
to be too small or too large. This modification has been found to improve the con- 
vergence when the desired transfer time interval is near a "knee'' in the At vs. x 
curve. 
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4. 3 Conic State Extrapolation as a Function of 
Transfer Angle - Theta Routine 

As with the Kepler Routine, the universal formulation of Battin in terms of 
the universal eccentric anomaly difference x is used in the Theta Routine. How- 
ever, in this case the functions a) may be expressed directly in terms of 

the specified transfer angle. The relationships used in this formulation are 

Uo(-f-;a) f ^ 

Tj^U;a) sin^-f- 

However, to evaluate these functions requires the final position magnitude, but 
from the position magnitude extrapolation equation it can be shown that 

2(cos I ( cos|-;^sin|- ) +-^ sin| ] -1 

Thus, from these quantities the extrapolated state vectors are calculated using the 
formulae in Section 4. Furthermore, these quantities may also be used to evaluate 
the transfer time interval from the generalized Kepler equation. 
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5 


DETAILED FLOW DIAGRAMS 


This section contains detailed flow diagrams of two Conic State Extrapolation 
Routines (Kepler and Theta) and the subroutines used by them. 

Each input and output variable in the routine and subroutine call statements 
can be followed by a symbol in brackets. This symbol identifies the notation for the 
corresponding variable in the detailed description and flow diagrams of the called 
routine. When identical notation is used, the bracket symbol is omitted. 

The Kepler routine utilizes the following subroutines, which are diagrammed 
in Section 5. 3; 

• Kepler Transfer Time Interval (Section 5, 3, 1) 

• Uj Series Summation (Section 5. 3. 2) 

• Q Continued Fraction (Section 5. 3. 3) 

• Kepler Iteration Loop (Section 5. 3. 4) ^ 

• Secant Iterator (Section 5. 3. 5) 

• Kepler Transfer Time Interval (Section 5. 3. 1) 

♦ Uj Series Summation (Section 5. 3. 2) 

• Q Continued Fraction (Section 5. 3, 3) 

• Extrapolated State Vector (Section 5. 3. 6) 

The Theta routine utilizes the following subroutines, which are diagrammed 
in Section 5. 3; 

• Q Continued Fraction (Section 5. 3. 3) 

• Extrapolated State Vector (Section 5. 3. 6) 
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UNIVERSAL PROGRAM INPUT 

CONSTANTS CONSTANTS VARIABLES 



Figure 3a. Kepler Routine, Detailed Flow Diagram 
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Figtire 3b. Kepler Routine, Detailed Flow Diagram 
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.2, 


^ < X ■‘CX 

min guess ma^ 


+ X 

^ min max 


guess 


Call Kepler Transfer Time 
Interval (Section 5^ 3^ 1) 

Input: X s[x^j.g].Jo*“ 

^ R k 

^ w ' max ' max 

0»<P“«- 'i'Vess I^^Wl 


-At <Al, 


^min "^^last ^gues^^^^ 
'^^.ast'^ ^*gties s 


“"gu^s last ^ ma 
^^guess ^ ^^ast 
max 
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Figure 3d. Kepler Routine, Detailed Flow Diagram 
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UNIVERSAL PROGRAM INPUT 

CONSTANT CONSTANTS VARIABLES 



Figure 4a, Theta Routine, Detailed Flow Diagram 
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Figure 4b. Theta Routine, Detailed Flow Diagram 
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V 



Figtire 4c. Theta Routine, Detailed Flow Diagram 
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5. 3 Subroutines Used by the Transfer Time and/ or the 
Transfer Angle Conic Extrapolation Routines 


5, 3, 1 Kepler Transfer Time Interval Subroutine 


INPUT VARIABLES 



Figure 5. Kepler Transfer Time Interval, Detailed Flow Diagram 
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Figure 6. Uj^ Series Summation, Detailed Flow Diagram 
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5. 3. 3 Q Continued Fraction Subroutine 


This subroutine is identical to the one used in the Conic Required Velocity 
routines. It should be noted that this routine has been designed to give double 
precision (16 digits) accuracy. 

INPUT VARIABLES 



Figure 7a. Q Continued Fraction, Detailed Flow Diagram 
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The symbol [x j denotes the largest integer whose magnitude 
does not " exceed the magnitude of x^ . 


Figure 7b. Q Continued Fraction, Detailed Flow Diagram 
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5, 3. 4 Kepler Iteration Loop Subroutine 

This subroutine is similar, though not identical, to the iteration loops used 
in the Conic Required Velocity routines. All could easily be combined into one 
routine but have been diagrammed spearately for purposes of clarity. 


INPUT VARIABLES 



Figure 8a. Kepler Iteration Loop, Detailed Flow Diagram 
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Call Kepler Transfer Time 
Interval (Sec. 5. 3. 1) 

Input: [x 

^ guess ^ arg^ 0 

€ , 3 , k 

w max max 



Figure 8b. Kepler Iteration Loop, Detailed Flow Diagram 
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5* 3* 5 Secant Iterator Subroutine 


routines. 


This subroutine is identical to the one used in the Conic Required Velocity 

. 

INPUT VARIABLES 

^t* ^ error' ^guess' '^^guess' 

*min’ ^%nin* *max' 


At ' . 
min 

At' 

max 


A'r - A? . 

guess min 

At -At 

guess max 


X or 

x|At» < e! 

' mav ' 


Ax = 0 


<0r 

V error ^ 


X " X \ 

guess minK g^ 

At' . / ^error 

min ' 


= lAt 


Yes >^x « -iN 

guess 

< rv ^ /%/ 


'*gue^s 

V Ax > X . / 
\ min^ 


( r^ rv V 

Vess ' ^maxU 

^^max r error 


X "X 

guess min 

At ' . 
min 


^min ” ^guess 

A^ . = ^ 
mm guess 


X ^ 

max 

^'^max 


X 

guess 

guess 


OUTPUT VARIABLES 


Ax, x^.^. At . , X ^ , At 

min' min max max 


Figure 9. Secant Iterator, Detailed Flow Diagram 
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5. 3. 6 Extrapolated State Vector Subroutine 


INPUT VARIABLES 



Figure 10. Extrapolated State Vector, Detailed Flow Diagram 







6. SUPPLEMENTARY INFORMATION 

The anal 3 ?tic expressions for the generalized Kepler equation and the extra- 
polated position and velocity vectors are given in Battin (Ref. 3 ). Battin also de- 
vised the normalized formulation and organized the overall flow diagram structure. 
However, the details of the structure and the testing were done by Shepperd and 
Robertson. 

A Newton iteration scheme may be devised at some time in the future. This 
would result in faster convergence near the solution point, but its worth depends 
on the expense of the derivative evaluation. It may be shown that if the derivative 
evaluation by itself takes more than 44% of the computation time used for the other 
calculations in one pass through the loop, then it is more efficient timewise to use 
the secant method. 
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